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We study the propagation of non-interacting polariton wavepackets. We show how two quali¬ 
tatively different concepts of mass that arise from the peculiar polariton dispersion lead to a new 
type of particle-like object from non-interacting fields—much like self-accelerating beams—shaped 
by the Rabi coupling out of Gaussian initial states. A divergence and change of sign of the diffusive 
mass results in a “mass wall” on which polariton wavepackets bounce back. Together with the Rabi 
dynamics, this yield propagation of ultrafast subpackets and ordering of a spacetime crystal. 


Field theory unifies the concepts of waves and par¬ 
ticles [1]. In quantum physics, this brought at rest the 
dispute of the pre-second-quantization era, on the nature 
of the wavefunction. As one highlight of this conundrum, 
the coherent state emerged as an attempt by Schrodinger 
to prove Heisenberg that his equation is suitable to de¬ 
scribe particles since some solutions exist that remain lo¬ 
calized [2]. However, the reliance on an external potential 
and the lack of other particle properties—like resilience 
to collisions—makes this qualification a moot point and 
quantum particles are now understood as excitations of 
the field. The deep connection between fields and par¬ 
ticles is not exclusively quantum and classical fields also 
provide a robust notion of particles, most famously with 
solitons [3]. The particle cohesion is here assured self- 
consistently by the interactions, allowing free propaga¬ 
tion and surviving collisions with other solitons (possi¬ 
bly with a phase shift). For a long time, this has been 
the major example of how to define a particle out of a 
classical field, until Berry and Balazs discovered the first 
case of a similar behaviour in a non-interacting context: 
the Airy beams [4] . These solutions to Schrodinger equa¬ 
tion (or equivalently through the Eikonal approximation, 
to Maxwell equations) retain their shape as they propa¬ 
gate as a train of peaks (or sub-packets) and also exhibit 
self-healing after passing through an obstacle [5]. The 
ingredient powering these particle behaviours is phase¬ 
shaping, assuring the cohesion by the acceleration of the 
sub-packets inside the mother packet. The solution was 
first regarded as a mathematical curiosity as it is not 
normalizable, till a truncated version was experimentally 
realized and shown to exhibit this dramatic phenomenol¬ 
ogy but for a finite time [6] . The Airy beam is now a rec¬ 
ognized particle-like object, in some cases emerging from 
fields that quantize elementary particles [ 7 ], thus behav¬ 
ing like a meta-particle. It is in fact but one example of 
a full family of so-called “accelerating beams” [8], that 
all similarly endow linear fields with particle properties: 
shape-preservation and resilience to collisions. 

In this Letter, we add another member to the family of 
mechanisms that provide non-interacting fields with par¬ 
ticle properties. Namely, we show that two coupled fields 
of different masses can support self-interfering wavepack- 



FIG. 1: (Golor online) a) Polariton dispersions. In red: the 
parabolic dispersion of the cavity photon, and the bare exci- 
ton. In blue: the polariton branches E± . The vertical dashed 
lines at A and Z 2 mark the inflexion points of the LPB. b) 
Effective masses for the LPB as a function of momentum: in 
purple the inertial mass mi, in blue the diffusive mass m 2 
(negative when A < /c < 22 ), in green the group velocity v-. 


ets, resulting in the propagation of a train of sub-packets, 
much like the Airy beam, but without acceleration, fully- 
normalizable and self-created out of a Gaussian initial 
state. Such coupled fields can be conveniently provided 
in the laboratory by polaritons [9] , the quantum superpo¬ 
sition of the spatially extended light and matter 

fields '0x(^,t), cf. Fig0 They find their most versatile 
and tunable implementation in semiconductors where ex- 
citons (electron-hole pairs) of a quantum well are coupled 
to the photons of a single-mode of a microcavity. We will 
consider the simplest ID case, realized for instance in 
quantum wires m but similar results hold for the more 
common planar geometry. Since polaritons can form con¬ 
densates giving rise to a wavefunction that describes their 
collective dynamics m, they are a dream laboratory to 
investigate the wavepacket propagation in a variety of 
contexts m, such as propagation of spin m, bullets m 
or Rabi oscillations m with technological applications 
already in sight [161 HZ]- Polaritons are highly valued 
for their nonlinear properties due to the particles self¬ 
interactions m, illustrated by a whole family of solitons 
(bright, dark, composite...) [TM22] . Recently, however, 
also the non-interacting regime has proved to be topical, 
with reports of skyrmions analogues [23] , band structure 
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engineering [24] focusing and conical polariton diffrac¬ 
tion [25|, internal Bosonic Josephson junctions [26|, em¬ 
ulates of oblique dark and half solitons [27] , Z topological 
insulator [28] or the implementation of Hebbian learning 
in neural networks [29] to name a few but illustrative 
examples. In most of these cases, interactions bring the 
physics to even farther extents rather than spoiling the 
underlying linear effect, that remains nevertheless the 
one capturing the phenomenon. The linear regime can 
be achieved at low densities m since the polariton in¬ 
teraction at the few particles level is small. In this case, 
the dynamics of the wavefunction l^ip) is ruled by the po¬ 
lariton propagator 11 such that |'0(t)) = Il{t — to) |'0(to)). 
In free space, the propagator is diagonal in k space m- 


{k'\Il{t)\k) = exp 


I ^ 


5{k-k') ( 1 ) 


where me is the photon mass, mx the exciton 
mass, A their detuning and their Rabi cou¬ 

pling. The eigenstates of the propagator, 11 (t) \\k))^ = 
ex.p{—iE±t) ||A^))j_, define both the polariton disper¬ 
sion E± = hk^mj^ -h 2A ^ /cq and the canonical polari¬ 
ton basis \\k))^ oc {E±{k),l)'^\k) where m± = {me ± 
^x)/(^c^x) are the reduced relative masses, kQ = 

^h?k^m^_ — 4:hk‘^Am- E 4(A^ -1- 40^) the dressed mo¬ 
mentum and \k) the plane wave of well-defined momen¬ 
tum k. We use the notation ||))^ for upper (+) and 
lower (—) polaritons. A general polariton state is thus 
expressed as \\ip))^ = f^(/>±(k) ||/^^))^ dk where 0±(/c) 
is the scalar-field polariton wavefunction. Except for a 
well-defined polariton state in /c-space, i.e., a fully de¬ 
localized polariton in real space, the photon and ex¬ 
citon components of a polariton cannot be jointly de¬ 
fined according to a given wavepacket 0(/c). Indeed, ex¬ 
cept if 0(/c) = S{k), one component gets modulated by 
the E±{k) factor needed to maintain the particle on its 
own branch. One striking consequence of this compos¬ 
ite structure is that a polariton cannot be localized in 
real-space, in the sense that both its photon and exciton 
components be simultaneously localized. Choosing (l){k) 
such that either 2 pe{x,t = 0) or = 0) is 6{x) re¬ 

sults in smearing out the other component in a pointed 
wavefunction surrounding the singularity of the localized 
field, as shown in Fig. |^a-b). Such constrains result 
in a rich phenomenology when involving a large enough 
set of momenta which, to the best of our knowledge, re¬ 
mained up to now safely hidden behind the simplicity 
of the problem. We devote the rest of the text to some 
of these remarkable effects, arising from the self-shaping 
and self-interferences of polaritons due to their composite 
structure, always in a non-interacting context. 

It has long been known that the mass imbalance me 
mx gives rise to a peculiar dispersion relation for the 
upper (E+) and lower (E_) polariton branches, shown 
in Fig. I^a) in blue along with, in red, the parabolic 
dispersions of the light photon and the heavy exciton. 



FIG. 2: (a) Localizing a polariton in space is possible for 
one of its component only (here the one in dashed red); the 
other field smears out to keep the particle on its branch, (b) 
Counterpart of (a) in energy-momentum space with forced 
localization of the photon field by the delocalized exciton. 
(c) Spacetime evolution of \'ipc{x,t)\‘^ and |'0x(ic,t)l^ with as 
an initial condition a photon of momentum ko = 0.5 /pm. 
The Bloch Sphere shows the quantum state trajectories along 
the line of the density plot in (c). (d) Configuration with 

A = — preventing splitting of the beam and resulting in 
ultrafast, Rabi-powered, propagating sub-packets, as shown 
for three snapshots of time in (e). 


meeting at /c = 0 (A = 0). The dynamics of a Gaus¬ 
sian wavepacket that is large enough in space to probe 
only parabolic portions of the dispersion in reciprocal 
space is essentially that expected from Schrodinger dy¬ 
namics m, diffusing with mean standard deviation of 
the packet size [32] : 

cTxit) = \/o-s(0) + (?it/[2m2CTj;(0)])2 . (2) 

Exciting one field only rather than eigenstate- 
superpositions (the polaritons) yield Rabi oscillations. 
Even in this simple case, there are subtleties brought 
by these non-parabolic dispersions. In particular, the 
degeneracy is lifted for some of the various concepts of 
masses, famously unified for the gravitational and inertial 
masses by Einstein as part of his theory of gravitation. 
For wavepackets, there are two different effective masses 
mi and m 2 [33], describing respectively propagation and 
diffusion. A wavepacket propagates with a group velocity 
= dkE±{k). This defines the inertial mass mi that 
determines the wavepacket velocity from de Broglie’s re¬ 
lation p = hk and the classical momentum p = mv± as: 


mi{E,k) = h^k{dkE)-'^. 


( 3 ) 
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FIG. 3: (Color online) (a-d) Propagation of lower polariton packets for various momenta and size, showing the emergence 
of the SIP for narrower packets, (e) Dynamics of an even narrower packet (cra;(0) = 2 pm with no momentum /cq = 0). (f) 
Current probability j at early times, showing the coexistence and interleaving of net counter-propagating flows, (g) phase map 
in a selected region, showing tt jumps associated to each sub-peak, (h) intensity profile at t = 90 ps with the evolution of the 
quantum state on the Bloch sphere corresponding to the path (from green to red) plotted in (e). Wavelet decomposition of (e) 
at t = 100 ps, and (j) in the same configuration but exciting with a momentum ko —12- (k) Spacetime honeycomb lattice when 
combining the SIP with Rabi oscillations by starting with a photon as an initial condition. (1) Zoom of the hexagonal lattice. 


A second mass m 2 , that we will call it the diffusive mass^ 
is associated with the spreading of the wavepacket ac¬ 
cording to Eq. (© and depends on the branch’s curva¬ 
ture; it reads: 

m2{E,k) = h\dlE)-K (4) 

These two masses mi, m 2 and the packet velocity V- 
are plotted in Fig. [^b) for the Lower Polariton Branch 
(LPB). Unlike parabolic dispersions, where they are 
equal, polariton dispersions yield qualitatively differ¬ 
ing inertial and diffusive masses. In particular, the k- 
dependent inertial mass mi imposes a maximum speed 
for the lower polaritons [32]. Beyond the inflexion 
point U, polaritons slow down if one increases their mo¬ 
mentum (at very large k the polariton becomes effec¬ 
tively a bare particle again with no such kinematic re¬ 
striction). The coupling of the two fields with differ¬ 
ent masses results in the heavier one lagging behind 
the other, as seen in Fig. |^c) for the case where a 
Gaussian photon wavepacket is imparted with a momen¬ 
tum k = 0.5/h achieved experimentally by send¬ 

ing a pulse at an angle and overlapping both branches. 
This prevents the photon and exciton packets to prop¬ 
agate Rabi-oscillating, and instead force a splitting in 
two beams—the orthogonal polariton states which are 
eigenstates for the corresponding wavevector, as shown 
by their trajectory on the Bloch sphere—connected by a 
Rabi oscillating tunnel. The Rabi oscillations only take 
place when there is a spatial overlap between the polari¬ 
tons. The two propagating packets maintain their co¬ 
herence despite their space separation and would Rabi 
oscillate if meeting again, due, for instance, to a ping- 
pong reflection [34]. The splitting in two beams can 
be minimized by tuning parameters to equalize the po¬ 
laritons masses, in particular the inertial ones. Com¬ 
bined with the bending of the Rabi oscillations in space¬ 


time, which can be achieved at nonzero detuning, this 
leads to propagation of Rabi oscillations, that produce 
ultrafast subpackets moving inside a mother packet, as 
shown in Fig.j^d) and for three snapshots in time in (e). 
The subpackets, continuously formed in the tail of the 
mother packet, propagate inside one order of magnitude 
faster, powered by Rabi oscillations, before dying in the 
head. Each sub-peak acquires properties of an identifi¬ 
able object that can be tracked in time. The full dynam¬ 
ics is available in an accompanying video [32]. Now on 
the diffusive mass m 2 : it diverges at the two inflexion 
points ii ^2 of the LPB and is negative in between. Ex¬ 
citing at the inflexion points thus cancels diffusion of the 
wavepacket as seen in Eq. Q and in Fig. [^a,b) with the 
propagation of a broad (crj(0) = 20 jam) lower-polariton 
wavepacket with an imparted momentum of (a) /cq = 0 
and (b) ko = ii. The excitation around the inflexion 
point has already been used to generate bright solitons 
and soliton trains [m [201 Ea ESI- In these cases, the 
soliton mechanism is the conventional interplay between 
negative effective mass and repulsive nonlinear interac¬ 
tions. The role of the high effective mass close to the 
inflexion point, which cancels the diffusion, was not how¬ 
ever fully estimated. 

The interesting phenomenology discussed so far illus¬ 
trates isolated features of the polariton propagation. A 
new physical picture emerges when combining several as¬ 
pects within the same wavefunction, leading to the con¬ 
cept of self-interfering packets (SIP). This is obtained 
when reducing the size of the wavepacket in real space, 
that is, increasing the staggering on the dispersion in 
momentum-space (cr/c(0) = l/cra,(0)), to an extent enough 
to probe polaritonic deviation from the parabolic disper¬ 
sion. In this case, the negative mass plays an explicit 
role. Negative masses are a recurrent theme in physics 
but this is typically meant for the inertial mass [37] . The 
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sign of the diffusive mass would seem, at first, not to 
play a role since it enters as a square in Eq. Q, and 
this is indeed the case for a regular packet with mo¬ 
mentum ii < k < i 2 . When straddling over the di¬ 
vergence, however, self-interferences occur between har¬ 
monics of the packet subject to the positive mass and 
others to the negative mass. This result in a complete 
reshaping of the wavefunction, as shown in Fig. gc,d) 
decreasing crx{0) down to 10 pm and 4 pm. The part of 
the packet that goes beyond the divergence is reflected 
back and interferes with the rest of the packet that still 
propatages forward, resulting in ripples. Reducing the 
packet to a^O) = 2 pm produces the striking pattern 
seen in Fig. [^e), without even the need of an imparted 
momentum. While for a parabolic dispersion, squeezing 
the packet in space merely causes a faster diffusion, in 
the polariton case, there is thus a critical diffusion be¬ 
yond which the packet stops expanding and folds back 
onto itself. Since this happens when the wavefunction 
encounters the inflection point of the polariton disper¬ 
sion, there is a “mass wall” against which the packet 
bounces back. If the dispersion also features another in¬ 
flection point at larger k, which can be the case for small 
enough exciton masses, this reflection happens again, this 
time resulting in a shielding from this self-interference 
of the core of the mother packet, as shown on the cut 
in intensity Fig. ih) (the diffusion cones are the so¬ 
lution of dlE_ = 0, cf. Fig. [^e)). More importantly 
from a conceptual point of view, as a result of this co¬ 
existence of masses of opposite signs within the same 
packet, the mother wavepacket {ip) fragments itself into 
two trains of daughter shape-preserving subpackets which 
travel in opposite directions. The overall momentum 
('0|p|'0) = 0 is null but the self-shaping of the wave- 
function redistributes it through its subpackets as a se¬ 
ries of nonzero momenta. Each sub-peak can be identi¬ 
fied by as a polariton as seen by following its quantum 
state on the Bloch sphere, that lies onto the meridian 
between |L) and |X), as shown in Fig. [^f) following a 
path (at t = 50 ps) from the central aera—shielded from 
the self-interferences—to the edge of the packet. The 
SIP can therefore be seen as a train of successive po¬ 
lariton packets, “emitted” by the area shielded from the 
self-interference at the rate of Rabi oscillations, and that 
retain their individuality as they propagate inside the 
mother packet. The full quantum state dynamics along 
these paths can also be seen vividly in a supplementary 
video [32]. Successive peaks furthermore feature a max¬ 
imal phase shift of tt in the phase 0(x, t) of the total 
wavefunction = |'0(x, t)| exp(i(/)(x, t)), as shown 

in Fig. [^g). Baring the fact that they do not involve 
self-interactions to account for their cohesion and other 
properties making them particles lookalike, these propa¬ 
gating subpackets behave in many respects as soliton-like 
objects. The analogy with Airy beams is conspicuous. 

One can gain additional insights into the nature of the 
SIP through the current probability j = ih/2mi{ip''dx'(p — 
xpdxtp*), in Fig. [sjf) where the packet is plainly 


seen to alternate backward and forward net flows. 
More sophisticatedly, considering the wavelet trans- 
form (WT) [SH] ^a,bW = (^) da;, 

in our case, of the Gabor wavelet family = 

^^exp(iu;x) exp(—x^/2), allows us to decompose the 
wavefunction into Gaussian packets, which are the ba¬ 
sic packets as far as propagation and diffusion are con¬ 
cerned. Such an extension of the Fourier transform is 
common in signal processing but has found so far little 
echo to study the dynamics of wavepackets [39] . We show 
in Fig.^i) the energy density |W/c,a:P of the wavefunction 
in the [x — k) plane at t = lOOps. One can see clearly how 
the self-interferences force the polariton packet to remain 
within the diffusion cone (blue dashed lines), by divert¬ 
ing the flow backward, (i) one or (j) two times when the 
second inflexion point {ko = ^ 2 ) is reached. Other funda¬ 
mental connections could be established. For instance, 
patterns strikingly similar to Fig. i: e)) were observed 
in the quenched dynamics of a quantum spin chains with 
magnons [40], a completely different system. This sug¬ 
gests that coupled light-fields feature fundamental and 
universal dynamical evolutions. Gombining this charac¬ 
teristic pattern with that of Rabi oscillations leads to the 
space-time propagation presented in Fig. [^c). The pro¬ 
tected area simply exhibits Rabi oscillations. The outer 
area is propagating upper polaritons and is not affected 
either by oscillations nor interferences. In the SIP area, 
however, sitting between the two mass walls, the inter¬ 
play of Rabi oscillations and self-interferences produces 
and hexagonal lattice. Here, instead of the emergence of 
propagating particles, a spacetime crystal is formed with 
the manifest ordering of the previously freely propagat¬ 
ing train of polaritons. This striking structure is, again, 
sculpted self-consistently out of a simple Gaussian state 
by the dynamics of coupled non-interacting fields. 

In conclusion, we have shown the intricate wavepacket 
propagations of coupled fields (polaritons). While the 
boundless diffusion of a Schrodinger wavepacket in a 
parabolic dispersion ultimately leads to complete in¬ 
determinacy, the polariton case can sustain traceable 
objects with always well-defined properties, such as 
their shape, position, momentum and quantum state. 
This gives rise to a concept of particles similar to 
that brought by the soliton in nonlinear media or Airy 
beams in non-interacting ones. While these are formed 
by self-interaction and phase-shaping, the individuality 
of polaritons is acquired and maintained through self¬ 
interferences powered by the Rabi coupling. This shows 
that even in the linear regime, the polariton dynamics 
is rich and able to produce intricate structures out of 
mere Gaussian initial states. This could lead to appli¬ 
cations, by imparting momentum powered by the Rabi 
oscillations or in the limit of few particles, for quantum 
computing, by a proper wiring and directing of the sub¬ 
packets, since all this happens in a strict linear regime. 
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Self-interfering wavepackets 
Supplementary Material 


INTRODUCTION AND OUTLOOK 


The theoretical model describing self-interacting 
wavepackets (SIP) is simple: two coupled ID Schrodinger 
equations in the linear regime (see section IV for the 2D 
case). The simplicity of the equation is no guarantee 
that its extent and depth are rapidly exhausted, as illus¬ 
trated by the Schrodinger equation, one of the most fun¬ 
damental equation of modern physics, for which the self- 
accelerating solution was discovered only in 1979. There 
should be hope, however, that some closed-form solutions 
are available. In Section III we show that the SIP, if it 


has such analytical expressions, does not seem to be re¬ 
ducible to a simple closed form. We could express it as a 
complex combinatorial superposition of Bessel functions, 
as could be expected from propagating packets, weighted 
by polaritonic factors, such as the dressed momentum k^. 
Before discussing this structure, we start in Section [ll| 
with more details on the exact solutions, both from the 
formalism and numerical simulations point of view, con¬ 
trasting in particular the multitude of ways that one has 
to poke at the polariton wavefunction. Last Section, |V| 
gives an overview of the rest of the Supplementary ma¬ 
terial, that consists of animated movies that illustrate, 
maybe better than equations, the mesmerizing dynamics 
of polaritons. 



FIG. 4: (Color online) Polariton propagation of delocalized 
wavepackets, as seen through the photonic ("ipc) and exci- 
tonic ('0x) components for the cases of (a) upper row: a 
lower polariton at t = 0 and (b) lower row: a photon. The 
diffusion is negligible over the selected time window and the 
staggering on the dispersion too small to evidence polariton 
self-interference effects. Parameters: Qr = 2 meV, me = 
0.025 meV ps^ , mx = 0.2 meV ps^ , cr^ = 20 pm. 


II. MORE POLARITON PROPAGATION 

The polariton propagator, Eq. (1) of the main text is 
easily found in k space (we work here with A = 0): 


{k'\n\k) = 


exp( 


4 


)(■ 


i40R 


/ k(^t \ I 

cos(^) -h I 




|p-m_ sin(-^ 


sp( 


i40R 

\ ‘ / k jt \ / — %k TTL-i-t \ f //Cq^\ ‘ k"^ * 

) sm{-^) exp(-pos(-g-) - sm{-f- 


exp( 


—ik^rrij^t 

4 


)sm(^) 



where we remind the important dressed momentum vari¬ 
able k^: 

kQ = ^k^m^ + , (S2) 

with m± = (me ± ^x)/(^c^x) (note again that k^ is 
a function of k). We also assume h = 1. Polaritons are 
maybe best formally defined as the states with a well- 
defined momentum and, consequently, also energy. They 
satisfy: 

n(t) \\k))^ = exp{-iE±t) ||fc))j^ , (S3) 

and as such are expressed as: 

P))± « 1^). (S4) 

for + (resp. —) the upper (resp. lower) polariton, with E± 
the pivotal polariton quantity, the dispersion: 

E± = k^rrij^ T • (S5) 


States (S4) form a canonical basis out of which a gen¬ 
eral polariton state is obtained by linear combination: 


/ oo 

<t^±m\k))±dk, (S6) 

-C» 

with 0±(/c) the scalar-field (upper/lower) polariton wave- 
function. All the results in this text follow from the 
impossibility to evolve such a general polariton state in 
time with the complex rotation of free propagation as 
in Eq. ( |S3| ), due to its two-component character. We 
have already discussed in the main text some important 
constrains of such a structure: for instance, that except 
for a well-defined polariton state in /c-space, i.e., a com¬ 
pletely delocalized polariton in real space, the photon and 
exciton components of a polariton cannot be jointly de¬ 
fined according to a given wavepacket 0(/c), e.g., a Gaus¬ 
sian packet, since one component gets modulated by the 
E± factor needed to maintain the particle on its branch. 
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FIG. 5: (Color online) Propagation of sharp packets. As seen in both real-space as a function of time (upper two rows) and 
through their dispersion (energy as a function of momentum, lower two rows) through the light emitted by the cavity (C) or 
through the direct exciton emission (X) for the cases of a lower polariton propagation (two columns on the left) and of a bare 
state propagation (two columns on the rigth). In the case of polaritons, (a-d), the exciton (upper row) or the photon (lower row) 
is perfectly localized at t = 0 with the other component defined such that the the particle remains on its branch. In the case of 
bare particles (e-h), the initial state is simply a photon (upper row) or an exciton (lower row) with the other component empty 
at t = 0. Greek-numbered panels correspond to the latin-numbered ones. Parameters: Qr = 0.5meV,mc = 0.5 meV ps^ 


Gaussian packets for both the photon and the exciton 
result in populating both polariton branches. The gen¬ 
eral case obviously admixes the two types of polaritons: 
l^/i) = f 0cr(^) ||^))± dk. These results that im¬ 

pose strong constrains on a polariton wavepacket must 
be contrasted with the conventional picture one has of 
the polariton as a particle, which is that of states \\k))^ 
and is, in good approximation, recovered for large enough 
packets as shown in Fig. The particle is here broad 
enough in space to have a small diffusion, cf. Eq. (2), 
and Rabi oscillations may be present depending on the 
state preparation, which however do not result in qualita¬ 
tive novelties. The situation is completely different when 
considering overlaps in k space, that is, sharp packets 
in real spaces. The propagation of such sharp packets 
is shown in Fig. for various initial conditions: as a 


lower polariton (left part of the figure) or as a bare par¬ 
ticle (exciton or photon, right part of the figure) and 
such that the excitonic component is perfiectly localized 
at t = 0 (upper row) or the photon component is (lower 
row). The wavepacket evolution is also shown as seen ei¬ 
ther through its photon or exciton field. Experimentally, 
the photon field is typically the one observed by record¬ 
ing the light leaked by the cavity. One sees variations 
around the themes exposed in the main text, with more 
or less pronounced features in some of the configurations. 
The clearest effects are obtained for confined excitons and 
observation through the cavity is always a good vantage 
point. Eor each case of propagation in space-time, we also 
show in Eig. the corresponding dispersion, which 

is the double Eourier transform. This shows how, in¬ 
deed, the lower polariton only populates its own branch. 
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FIG. 6: (Color online) Propagation in /c-space for non¬ 
diffusing packets of, left column, a lower polariton and, 
right column, a photon, as seen through, from up to bot¬ 
tom: \ 2 pc{k,t)\‘^, and Parameters: 

= 0.5 meV, me = 0.5 meV ps^ . 


It also provides an alternative view of the various local¬ 
ized states, e.g., the localized photon does not populate 


the photon-like part of the polariton branch in the exci- 
ton spectrum (panel (5) and vice-versa with the localized 
exciton in the photon spectrum (panel a). 

We now proceed with the underlying Mathematical ex¬ 
pressions. For, say, the lower polariton case prepared so 
that the photon component is perfectly localized at t = 0 
(in the case where, for concision in the notation, we as¬ 
sume from now on mx ^ oc), we find from Eqs. (S1-S6): 


/ 7^2 _ \ 7^2 _ iL2 

(S7a) 

V’x(fc,i) = exp , (S7b) 


The result is extremely simple in this picture (time- 
momentum). There is no time dynamics for the den¬ 
sity |'0x,cP- the perfectly localized exciton at t = 0 
results in a completely delocalized wavefunction at all 
times IV^xP = 1- The corresponding photon wavefunc¬ 
tion is qualitatively similar to a Voigt lineshape in k- 
space, which we will use in next Section to derive ap¬ 
proximated expressions. It is shown on the left column 
of Fig. In all cases, however, there is of course a dy¬ 
namics of the wavefunction itself, as seen through its real 
and imaginary parts on the figure. There is a slowing 
down of the oscillations with increasing \k\. The Fourier- 
transform in k of this pattern gives the spacetime prop¬ 
agation in Fig. [^a). 

The case of only the perfectly localized exciton as 
the initial condition (i.e., in the photon vacuum rather 
than the field needed to provide a lower polariton), i.e., 
'0_(x,t = 0) = (J(x),0)^, is given directly by the 
columns of Eq. ([ST}: 







exp 



{—iftjit) sine 



(S8a) 

(S8b) 


The corresponding propagation in momentum-space is 
shown on the right column of Eig. There is, this 
time, a dynamic in the density, namely, the Rabi oscilla¬ 
tions, with, in contrast to dynamics of the real an imagi¬ 
nary parts, a speeding up of the oscillations with increas¬ 
ing |/c|, corresponding to the effective Rabi frequency of 
effectively detuned exciton-photon coupled states. The 
Eourier-transform in k of this pattern gives the space- 
time propagation in Eig. |^e). 

We have also discussed in the main text how a polari¬ 
ton wavepacket propagates with a group velocity v± = 


dkE±{k) that, for the lower polariton, features a local 
maximum. It reads: 

— — ^ 

k k 2m- 

V- = - , , (S9) 

2m+ r 4m_A , 4m^(40|^+A2) 

+ —p— 

which is shown on Eig. in blue. The local maximum 
is obtained at the first inflexion point of the dispersion, 
/c = ii, where the polaritons also do not diffuse. In¬ 
creasing the momentum makes the particle heavier and 
actually reduces its speed. A local minimum is attained 
at the second inflexion point, k = i 2 ^ where the polari- 
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FIG. 7: (Color online) Speed of the lower polariton as a func¬ 
tion of its imparted momentum: there is a maximum speed 
for lower polaritons. Larger momenta result in slower parti¬ 
cles. Local maxima are given by the inflexion points. At large 
enough k, the polariton becomes exciton-like and suffers no 
such restriction. If the exciton mass is infinite, polaritons 
have an absolute maximum velocity. 



FIG. 8: (Color online ) (a,b) |'0x(^c,t)r calculated fro m th e 
exact result Eq. (S7b) and the approximated one, Eq. (Sll). 
(c) Normalized intensity at the center of the wavepacket {x = 
0): exact result through the Eourier Transf orm o f Eq. dSTbt 
(blue line) and the analytical expression Eq S12 (red dashed 
line). The approximation appears to be exact at a: = 0. 


ton also propagate without diffusion but now with a low 
speed. At larger Eq. ( |S9[ ) becomes linear and tends 
to the speed of the bare exciton, as indeed for k ^ 0, 
V-{k) k/mx^ cf. Fig. in dashed red for a non de¬ 
tuned system. If there is no second inflexion point (for 
an infinitely heavy exciton mass), there is an absolute 
maximum speed for the polaritons since the bare exciton 
group velocity vanishes (green curve). 


III. APPROXIMATIONS 

We could not find manageable closed form expres¬ 
sions for the Fourier Transform of the term exp(/cQ — k‘^) 
that capture the (lower) polaritonic self-interference ef¬ 
fect. The term k‘^ — k^ for the (lower) polariton momen¬ 
tum distribution can however be well approximated by 
a Voigt distribution, since it combines both exponential 
and fat-tail types of decay. Since the fat-tail is expected 
to play a dominant role qualitatively, we assume simply 
a Lorentzian distribution f{k) = dOpmc/ (l + ^n^mc ) ’ 


thus approximating Eq. (S7b) by: 




1 + fc2/(4mcriR) 


(SIO) 


which Fourier Transform can be obtained by a series 
expansion of the exponential, providing the real-space 
dynamics of the SIP: 


tpx{x,t) ~ 

^ n!2 

n=0 

(Sll) 


where Ka{z) is the modified Bessel functions (solution 
of th equation z‘^y" + zy' — {z‘^ n‘^)y = 0) and with 

a = n + 1/2. The propagation of the wavepacket cal¬ 
culated numerically with the Fourier Transform (a) and 
the one obtain with the analytical formula in Eq. |S11| 
(b) are compared in Fig. There is a good qualitative 
agreement, with only the size of the propagation cone 
differing slightly, due to the width difference between 
the Lorentzian distribution and the actual one. One can 
therefore trust the approximation to give some insights 
into the nature of the SIP. First, the SIP is indeed a 
phenomenon of many interferences. The convergence of 
the series is obtained for a number of terms in the sum 
that increases linearly with time t, showing how each new 
peak arises from an added term and thus a next order in 
the interference. Also, at the center of the wavepacket 
{x = 0), the previous expression can be reduced to a 
simple form: 


^x{x = 0,t) cz -27r^flRmc e 2 ^ (S12) 


where Jn{z) is the Bessel function of the first kind. Since 
the departure between Eq. (Sll) and the exact numeri¬ 


cal solution is mainly due to the extent of the envelope 
of the momentum, one can expect a better agreement 
at X = 0, and indeed we find that there is a perfect 
match, as seen in Fig. [^c). Equation confirms 

that the successive peaks that shape the SIP appear at 
the Rabi frequency flp. The photon mass me (we remind 
that we assumed here an infinite exciton mass) acts only 
on the intensity. In the same way, one can obtain the 
corresponding series for the photon field: 


tpcix,t) ~ 

4^imcnR)^{imnr\x\^~^K0{2Vn^\x\) 


E 

n=0 


i!r(n) 


(S13) 


with /3 = n—1/2, showing how both structures are tightly 
related and the sort of complexity that describes them. 
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FIG. 9: (Color online) SIP in two-dimensions, (a) With 
full radial symmetry when self-interfering in all directions. 
A central disk is shielded from interference just like in the 
ID case and feature a flat plateau of lower polaritons. (b) 
When imparted with momentum, the packet self-interfere in 
the direction of its propagation and diffuses radially. Pa¬ 
rameters : Qr = 2meV,mc = 0.025meVps^mx = 
0.2 meVps^ /im“^, (a) cr^ = 2 pm, (b) cTx — ^ pm, kx — ky = 
ii. 


Series for photons or excitons as initial conditions can 
also be obtained from Eq. S8a and S8b[ involving Hyper¬ 
geometric functions which, however, are too cumbersome 
to be written here. 


IV. SIP IN TWO DIMENSIONS 

So far, we have considered ID cases, which are in¬ 
deed possible in heterostructures by confining in the two 
other dimensions (quantum wires). Polariton propaga¬ 
tion is however popular in the 2D geometry as well and 
we quickly discuss what happens in this configuration. 
Since the system is linear and uncoupled in transverse 
coordinates, the dynamics follows trivially from the pre¬ 
vious results and symmetry. With a full radial symmetry, 
all the phenomenology is conserved in 2D, with rotational 
invariance. This is shown in Fig. |^a) where the prop¬ 
agation of a SIP is shown after 100 ps, after preparing 
a narrow Gaussian wavepacket at t = 0 (and at the ori¬ 
gin of the plane). This is the counterpart of the case of 
Fig. 3(e) of the main text. The propagation and thus the 
shape of the wavepacket are similarly determined by the 
dispersion and the way it is excited. There is also, de¬ 
pending on the proximity of the second inflection point, 
a central area shielded from the interferences, which is 


now a disk, whose diameter is also determined by a zero 
of = 0. Also, like in the ID case, one can impulse 

the propagation of the packet in a desired direction by 
imparting a momentum. This is shown in Fig. [^b) with 
the propagation of a smaller wavepacket when exciting 
the dispersion at the first inflexion point = ky = ii). 
By using squeezed Gaussian packets, one can propagate 
a SIP with fronts that remain parallel and orthogonal to 
the direction of motion. 

V. MOVIES 


We also provide two movies that illustrate vividly the 
polariton dynamics. 

The movie I-QuantumState.avi animates Fig. 3 (e,h) 
of the main text. It shows how a narrowly squeezed lower 
polariton wavepacket self-interferes and produces, as a re¬ 
sult, two trains of sub-packets propagating back to back, 
emerging from a polariton sea shielded from the interfer¬ 
ences. To show how the overall structure of the SIP is 
connected to the individual identity of each packet, we 
also present dynamically on the Bloch sphere the evolu¬ 
tion of the quantum state on a path that links the center 
of the wavepacket to the side of the diffusion cone (from 
the green to the red point on the density plot, mached 
with the green and red arrows on the sphere). The state 
in the central area—shielded from interferences—is the 
lower polariton. In the interferences zone, crossing a 
fringe induces a loop on the sphere that crosses the merid¬ 
ian of states |C)-|L)-|X). Each peak converges in time 
towards a well defined polariton state on the meridian. 

The movie II-PolaritonRiffle.avi animates the 
case of Fig. 2(d,e) of the main text. It shows the propaga¬ 
tion of a photon wavepacket (at t = 0) with a momentum 
and negative detuning. A judicious choice of parameters 
permits to maintain a spatial overlap between the two 
polaritons, conserving the Rabi oscillations which, due 
to the detuning, are bent in spacetime (see main text). 
This results in the propagation of ultrafast subpackets 
within the mother packet. This is well seen on the ani¬ 
mation after 35 ps of animation time, which is the time 
needed for the packet to develop the structure (the initial 
condition is a Gaussian packet). 










